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We study graphene antidot lattices (GALs) in magnetic fields. Using a tight-binding model and 
a recursive Green's function technique that we extend to deal with periodic structures, we calculate 
Hofstadter butterflies of GALs. We compare the results to those obtained in a simpler gapped 
graphene model. A crucial difference emerges in the behaviour of the lowest Landau level, which 
in a gapped graphene model is independent of magnetic field. In stark contrast to this picture, we 
find that in GALs the band gap can be completely closed by applying a magnetic field. While our 
numerical simulations can only be performed on structures much smaller than can be experimentally 
realized, we find that the critical magnetic field for which the gap closes can be directly related to 
the ratio between the cyclotron radius and the neck width of the GAL. In this way, we obtain a 
simple scaling law for extrapolation of our results to more realistically sized structures and find 
resulting quenching magnetic fields that should be well within reach of experiments. 

PACS numbers: 81.05.uc, 78.20.Ls,78.20.Bh 



I. INTRODUCTION 

Semiconductor antidot lattices have revealed a range of 
intriguing transport phenomena. In particular, transport 
in magnetic fields has been applied in studies of magneto- 
resistance, localization, Shubnikov - de Haas oscillations, 
and quantum Hall efFectsjii^ Traditionally, such struc- 
tures have been based on GaAs heterojunctions. Re- 
cently, however, graphene antidot lattices (GALs) have 
been proposed^i^ and magneto-transport studies have 
been realized experimentally<^"— Whereas GaAs antidots 
are typically produced by dry etching or local ion implan- 
tation in molecular beam epitaxy grown heterostructures, 
GALs are produced by simply etching arrays of holes into 
graphene sheets placed on suitable insulating substrates. 
Etch masks can be fabricated using either e-beam lithog- 
raphy or block co-polymers. In this manner, features of 
the order 10-20 nm have been obtained^ii^. In addition 
to GALs fabricated via etching, the patterned hydrogen 
adsorption technique^i also produces structures resem- 
bling periodic antidots. Here, hydrogenated graphane- 
like islands form 'forbidden' regions similar to holes but 
on a few-nanometer scale. Very recently, direct e-beam 
writing of small GALs with holes around 2 nm has been 
demonstrated 

While magneto-transport is well documented in exper- 
imental GAL studies, the theoretical description of ex- 
tended GALs in magnetic fields is a challenging task, 
because the associated vector potential breaks trans- 
lational invariance. Structures comprising only a sin- 
gle or few antidots are more easily analysed. For in- 
stance, Yang and coworkers^^iii have included magnetic 
fields in studies of isolated graphene antidots and small 
arrays, with perforations modelled as circular electro- 
static potentials. The present authors studied an iso- 



lated antidot using a mass term barrierJ^ The prob- 
lem of modeling periodic arrays of graphene antidots in 
magnetic fields is much more involved, however. The 
energy spectrum of electrons in periodic potentials sub- 
ject to magnetic fields takes the form of self-similar Hof- 
stadter 'butterflies' Hofstadter butterflies have pre- 
viously been studied in pristine graphene ; ^^'^^ bilayer 
graphene twisted bilayer graphene^ graphene with 
point defecta^i and graphene quantum dotsi^ Very re- 
cently, Hofstadter spectra have been studied experimen- 
tally in single and bilayer graphene on hexagonal boron 
nitride<^i^ In all cases, an intriguing self-similar struc- 
ture emerges whenever the characteristic magnetic length 
becomes comparable to the geometrical period. The rela- 
tively large size of this period is challenging for an atom- 
istic description of the spectrum, however. In pristine 
graphene, a magnetic field leads to a gapped Landau 
level structure. Hence, for graphene modified in order 
to induce band gaps (such as GALs, electrically biased 
bilayerS)^ quantum dots?^ and geometrically sharp elec- 



trostatic gate! 



an interesting interplay between ge- 



ometric and _B-ficld induced band gaps is expected. A 
simplified model capturing the essentials of geometrically 
induced band gaps is that of gapped graphene^^ for which 
a unique Landau level structure^ and magneto-optical 
response^ have been predicted. For magnetic effects, 
though, little is known about the reliability of gapped 
graphene as an approximation to graphene with a gapped 
spectrum induced by geometric modifications. 

In the present work, we study the properties of GALs 
in magnetic fields using a tight-binding model with mag- 
netic effects included via a Peierls phase. We compare 
true antidot geometries to gapped graphene approxima- 
tions and find important differences. Most importantly, 
we observe band gap quenching by the magnetic field 



in GALs whereas in gapped graphene the band gap is 
completely robust for arbitrarily large fields. Hence, the 
behavior of GALs is reminiscent of graphene quantum 
dots2^ and graphene nanoribbonsi ^^''^^ We also demon- 
strate simple scaling laws to extrapolate our results to 
experimentally feasible structures. 
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A GAL is modelled as an infinite, periodic array of 
circular holes in the graphene sheet. The super lattice 
spanning the holes is assumed to have regular triangular 
synnetry, in order to ensure the existence of a full band 
gap for all unit cell geometries i^^iM We denote the radius 
of the hole and the side-length of the unit cell by R and 
L, respectively, both of which are measured in units of 
the graphene lattice constant a. Hence, a given GAL is 
designated by {L, i?}, where L (but not necessarily R) is 
an integer. 

We use a nearest-neighbor tight-binding model to 
study the energy spectrum of 7r-elcctrons within the 
GAL, i.e. the hopping integral tij is taken as —t with t = 
3 eV for neighboring sites and vanishes otherwise. The 
effect of the magnetic field is included via a Peierls phase 
added to the hopping term between atomic sites i and j : 
tij — >■ iye"^, with the phase given as = (e/h) A-dl. 
Here, and Rj denote the positions of the atomic sites, 
while A is the magnetic vector potential. The graphene 
sheet is in the xy plane and the magnetic field is taken 
to be constant and directed perpendicularly to the sheet. 
Using the Landau gauge A = Bxy the Peierls phase be- 
comes 
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{x,+Xj){yj - y,). 



(1) 



The crucial advantage of the Peierls phase approach is 
that lattice periodicity can be restored provided a suit- 
able " magnetic supercell" containing several original unit 
cells is constructed. To this end, we first construct a rect- 
angular unit cell containing two antidots. This doubled 
cell is then repeated N times in the x direction as shown 
in Fig. [T] Periodicity is subsequently ensured by requir- 
ing that the phase shift difference be an integer multiple 
of 2tt for a pair of neighbor sites at the left- and right- 
most ends of the supercell, i.e. separated by a distance 
of 3NLa. The smallest separation along the y axis for 
coupled atoms is and so the minimum B- field re- 

quired for periodicity is 
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V^eNLa^ ' 



(2) 



When the fiux $ = B^/Sa^ /2 through a graphene unit 
cell equals a flux quantum $o = h/e, the energy spec- 
trum has been restored to the unperturbed one. Writ- 
ing B — nB^in with n integer we find the relative fiux 
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FIG. 1: (a) Enlarged supercell used for calculating the mag- 
netic properties of GALs. The smallest rectangular GAL unit 
cell is repeated A*' times to ensure periodicity of the Peierls 
phase in the hopping terms of the tight-binding Hamiltonian. 
(b) Rectangular unit cells for two examples of GALs with 
similar neck widths. 



— We can thus calculate the Hofstadter spectra 
by varying n over the range [Q,NL\. In practice, ad- 
vantage is taken of the fact that the same spectrum is 
obtained for identical ratios n/N ^ and thus calculations 
can be performed on much smaller magnetic unit cells 
whenever n and TV are commensurate. 

From Eq. © it follows that the number of repeated 
unit cells in the supercell required to handle a cer- 
tain magnetic field is inversely proportional to the field 
strength. Because of the relatively large size of the fun- 
damental unit cell of a GAL, the number of atoms in the 
magnetic supercell becomes very large even for substan- 
tial magnetic fields. We thus cannot easily rely on stan- 
dard diagonalization techniques to determine the prop- 
erties of GALs in magnetic fields. Instead, we expand 
on well-known recursive Green's function methods.— 
These methods are commonly applied to transport stud- 
ies of finite or semi-infinite structures sandwiched be- 
tween semi-infinite leads. Here, periodicity of the central 
region is at most along a single direction and the Hamilto- 
nian can easily be written in block-diagonal form, which 
is the basis upon which the recursive Green's function 
formalism is built. In our case, periodicity is along both 
directions. We note that one immediate solution to this 
problem is to perform a transformation of the Hamilto- 
nian to force it into block-diagonal form, which is possi- 
bly even with periodic boundary conditions in both di- 
rections. This, however, requires a doubling of the size 
of the individual cells used in the recursive algorithm 
which, for large GALs, we find to be slower than our 
method for calculating the density of states. In the ap- 
pendix, we present our extension of the recursive Green's 
function formalism so that periodicity can be handled 
without reshuffling the elements of the Hamiltonian. We 
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FIG. 2: (Color online) Hofstadter butterfly for the gapped 
graphene model, with the gap set to match that of the {4, 2} 
antidot lattice. Shown is the density of states as function of 
energy and the relative flux through a graphene unit cell. To 
ease visibility, the logarithm of the density of states is plotted. 
The bottom panel shows a closer view of the low-flux part of 
the spectrum. The dashed lines show the lowest five Landau 
levels expected from a Dirac treatment of the problem. 



use this formalism to calculate the density of states ac- 
cording to p{E) = — 7r^^Im{Tr[G"']}, where the retarded 
Green's function — {E + ij ~ H)^^ includes a small 
broadening term 7. A related method based on stochas- 
tic evaluation of G"" has previously been used to study 
magnetic effects in tvifisted bilayersj^ 



III. HOFSTADTER BUTTERFLIES 

As mentioned above, a phenomenological gapped 
graphene model may be used as an approximation of 
the true GAL geometry. Here, alternating on-site en- 
ergies ±A are assigned to the atoms belonging to the 
two sublattices. This model produces a band gap of 
Eg = 2A, which is adjusted to the gap of the full GAL 
band structure. For the optical response, we have pre- 
viously demonstrated good agreement between the GAL 
and gapped models for photon energies close to the band 
gap<^ To compare with the magnetic spectra for GALs 
we show in Fig. [2] the Hofstadter butterfly of a gapped 
graphene model, wherein the effect of the antidot lattice 
is included via such a staggered potential. The figure 
shows the density of states as a function of energy and the 
relative flux through a graphene unit cell. As expected 
from the discussion above, we see that when ^/^o = 1 
the spectrum in absence of a magnetic field is restored. 
One striking characteristic that emerges in the figure are 
the zeroth Landau levels, which for gapped graphene 
sit at ±A, i.e. exactly separated by the band gapj^ 



These zeroth Landau levels are characteristic of graphene 
and exhibit energies that are independent of magnetic 
fieldj^"— Whereas for low magnetic fields, where the 
Dirac model is applicable, the zeroth Landau levels are 
indeed completely independent of magnetic field, we do 
see a broadening and eventual splitting of the levels into 
multiple bands as the magnetic field becomes sufficiently 
strong. This is equivalent to what is seen for ordinary, 
gapless graphene 4^ 

Note that the full range of relative magnetic fiux shown 
in the figure corresponds to unrealistically large magnetic 
fields up to 79 kT, well beyond the reach of experiments. 
A zoom of the low-flux part of the butterfly is shown 
in the bottom panel of Fig. [2l For comparison we also 
show the Landau levels expected from a Dirac equation 
treatment. Ens = syA^ + nfi^ajJ, with s = ±1 and n 
a non- negative integer. Here, ujc = V^v^/Ib, with vp 
the Fermi velocity of graphene and Ib = y^H/JeB) the 
magnetic length. Note how the fractal structure only 
really emerges once the relative flux is of the order of 
$/$o = 0.05, corresponding to a huge magnetic field 
strength of the order of B ~ 4 kT. The crucial param- 
eter for seeing signatures of the fractal structure is that 
the magnetic length should be of the order of the lattice 
constant of the material. In pristine graphene ~ and in- 
deed in any bulk material - this requires unrealistically 
large magnetic fields. Thus, the much larger superlat- 
tice introduced by the GAL is a way of overcoming this 
obstacle. A similar explanation applies to the recent ex- 
perimental studies of effects related to Hofstadter butter- 
fiies in single- and bilayer graphene on hexagonal boron 
nitride<^i^ Here, a moire superlattice is formed due to 
the slight mismatch in lattice constant between graphene 
and boron nitride. In the present work, the fundamental 
superlattice period is given by the GAL lattice constant 
\/3La. Thus, due to the B~^^^ behavior of Ib, the in- 
field, at which novel magnetic features become visible, is 
reduced roughly by a factor of L^. 

The Hofstadter butterflies for two examples of GALs, 
{4, 2} and {12, 8.6}, are illustrated in Figs.HandU The 
geometries of these GALs are illustrated in Fig.[Tl3. Note 
that the radius is chosen such that we ensure that no dan- 
gling bonds are created. While these structures, having 
feature sizes below one nanometer, are smaller than what 
can be achieved experimentally, we will demonstrate be- 
low that certain scaling laws can be used to extrapo- 
late results to GALs with more realistically sized fea- 
tures. Focusing first on the {4, 2} GAL we note that 
while the structure of the spectrum is significantly richer 
than that of gapped graphene in Fig. [21 many features 
of the gapped graphene spectrum are preserved. In gen- 
eral, the large regions devoid of eigenstates for gapped 
graphene tend to be connected by additional bands in 
the GAL case. This structure becomes increasingly rich 
as the size of the GAL unit cell is increased, as is evident 
in Fig. 01 Comparing the two GALs, one significant dif- 
ference is the additional states at zero magnetic field for 
the {12,8.6} GAL, marked by arrows in Fig. [H These 
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FIG. 3: Hofstadter butterfly for the {4, 2} GAL, showing the 
density of states as function of energy and the relative flux 
through a graphene unit cell. To ease visibility, the logarithm 
of the density of states is plotted. Note the quenching of the 
gap by the magnetic field around <I>/"l>o = 0.1, shown in more 
detail in the zoom in Fig. [5] 
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FIG. 4: Hofstadter butterfly for the {12, 8.6} GAL, showing 
the density of states as function of energy and the relative flux 
through a graphene unit cell. To ease visibility, the logarithm 
of the density of states is plotted. The arrows indicate the 
states emerging due to regions of zigzag edges at the edge of 
the antidot. 



emerge due to local regions of zigzag geometry at the edge 
of the hole, which tend to induce localized states4i The 
localized nature of the states residing in this band are re- 
flected in its behaviour as the magnetic field is increased, 
where compared to other energy bands it remains very 
narrow with little splitting of the energy levels. 
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FIG. 5: Closer view of the low-flux part of the Hofstadter 
butterflies for {4, 2}, {8, 5}, {9, 6} and {12, 8.6} GALs, show- 
ing the quenching of the gap by the magnetic fleld. Note that 
these GALs all have similar minimum neck widths. 



IV. MAGNETICALLY INDUCED BAND GAP 
QUENCHING 

The most striking difference between the results for the 
GALs in Figs. |3] and H] and those of the simpler gapped 
graphene model shown in Fig. [5] is the behaviour of the 
lowest Landau level as the magnetic field is increased. 
A Dirac treatment of gapped graphene predicts that the 
energy of this state should be Eg/2, independent of mag- 
netic field, a behaviour that is confirmed in Fig. [5] also 
for very large magnetic fields. However, the spectra for 
the GALs show an entirely different behaviour. Here, the 
band gap is quenched as the magnetic field strength is in- 
creased and eventually closes entirely. This is reminiscent 
of what is seen for armchair graphene nanoribbons^i^^ 
and graphene quantum dots.— The crucial difference be- 
tween GALs and a gapped graphene model is the addi- 
tional characteristic lengths introduced by the antidots. 
In particular, we propose that the gap in a GAL will 
be quenched once the magnetic length Ib is of the or- 
der of the minimum neck width of the GAL. In a simple 
picture, when the cyclotron radius becomes sufficiently 
small, the individual eigenstates do not sample the lattice 
sufficiently for the band gap to be fully resolved. Because 
of the unique property of the lowest Landau level in bulk 
graphene, namely that it sits at the Dirac point energy 
regardless of magnetic field strength, this will result in a 
diminished band gap. We stress that we do not expect a 
similar effect to occur in 2DEGs based on semiconductor 
hetcrostructures, where the energy of the lowest Landau 
level is proportional to the magnetic field strength. 

To illustrate the discussion above, we show in Fig. [5] 
zooms of the low-flux region of the spectra of different 
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FIG. 6; (Color online) The magnetic length Ib correspond- 
ing to the critical flux for which the band gap is closed as a 
function of the neck width ?i>/a ~ \/3i — 2_R of the GAL. The 
dashed lines shows a best fit of the form oIb = to. 



GALs, more clearly illustrating the magnetically induced 
band gap quenching. We show four examples of GALs, 
all of which have similar neck widths. Despite significant 
differences in the band gaps at zero magnetic field we find 
that the magnetic flux for which the gap is completely 
quenched is very similar for all four GALs. Note the ap- 
pearance of additional very narrow bands at zero field in 
the case of the {9, 6} GALs (near ±0.45 eV) and {12, 8.6} 
GALs (near ±0.25 eV). As discussed above, these states 
appear due to quasi-localized states residing predomi- 
nantly on regions of local zigzag geometry at the edge 
of the antidots. Ignoring these additional states for now, 
and taking the band gap to be defined as the gap between 
the wider energy bands, we see that the dependence of 
the band gap on magnetic field is very similar in all four 
cases. 

In Fig. |6] we illustrate the geometry dependence of 
the magnetic length Ib corresponding to the critical fiux 
$c, for which the band gap is completely closed, on 
the minimum neck width. The neck width of a given 
GAL is approximately w/a = \piL — 2R, an expres- 
sion that is obviously more accurate for larger struc- 
tures where the exact atomic details can be safely disre- 
garded. As expected from the physical picture described 
above, we find a clear linear dependence between the 
two quantities, such that the critical fiux is defined via 
the simple relation als = w, with a ~ 3.1. Note that 
Ib — 25.7 nm T^/^/ \/B, suggesting that the gap would be 
quenched at realistic magnetic field strengths for larger 
structures. Indeed, the scaling law provides us with a 
means of extrapolating our theoretical results to more 
realistically sized structures. As an example, we consider 
a GAL with a lattice constant of 60 nm and an antidot 
radius of 25 nm, which represent experimentally feasible 
feature sizes 4 The neck width of this GAL is approxi- 
mately If ~ 54 nm, from which the scaling law predicts 
a critical flux of — 2.82 x 10~'^<f>o. This translates to 
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FIG. 7: (Color online) Magnetic flux dependence for a {12, 5} 
GAL of (upper panel) the density of states (in arbitrary units) 
at the Dirac point energy and (lower panel) the overlap with 
the antidot edge of the eigenstate corresponding to the lowest 
(positive) energy. The overlap is shown relative to the overlap 
assuming uniform distribution of the eigenstate. See the main 
text for details. In both panels, the dashed line shows best 
flts of exponential growth or decay, respectively. The insets 
in the lower panel show the probability density for the two 
eigenstates marked with circles. 



a magnetic field strength of i? ~ 2.2 T, well within reach 
of experiments. 

To further illustrate the interplay between the two 
length scales set by the magnetic field and the neck width 
of the GAL, we have studied the eigenstate nearest to the 
Dirac point energy. In particular, we consider the over- 
lap J^^^^\ilj{x,y)\'^dxdy of the probability density with 
the edge of the antidot. In practice, we take this as be- 
ing equal to sum over the absolute square |c„p of the 
expansion coefficients of the 7r-orbitals of carbon atoms 
with only two nearest neighbours. The lower panel of 
Fig. [7] shows the edge overlap as a function of the rela- 
tive magnetic fiux for a {12,5} GAL. The edge overlap 
is shown relative to the overlap that would be expected 
if the probability density were evenly distributed across 
all carbon atoms, i.e. we plot I]neodgo hi]"^ i^tot / Ncdgc) , 
where A^tot is the total number of carbon atoms in the 
structure while A'edge is the number of carbon atoms with 
only two nearest neighbours. In the figure, we also in- 
clude a best fit to the data, showing a clear exponential 
decay of the edge overlap with magnetic flux. For com- 
parison, we show in the upper panel the density of states 
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right at the Dirac point energy versus the magnetic flux. 
Here, we see an exponential increase of the DOS as the 
eigenstate is pushed away from the edge of the antidot. 
The critical flux where the band gap is closed corresponds 
quite accurately to the point where the exponential de- 
cay of the edge overlap starts to level off (for larger mag- 
netic fields than included in the figure). In the insets 
of the lower panel of Fig. [3 we show two charge densi- 
ties corresponding to the magnetic flux and edge over- 
laps indicated with circles in the figure. Because of the 
high degeneracy of the eigenstates we show the probabil- 
ity densities summed over all states in the lowest band 
of energies. The black dots indicate carbon atoms while 
the size and shading of the filled, red circles illustrate the 
probability density. To ease visibility we show just the 
first four rectangular GAL unit cells, which correspond 
only to a small region of the full magnetic unit cell used 
in the calculations. The transition from a state with sig- 
nificant overlap with the antidot edge to a state localized 
predominantly between the antidots is quite evident. 

These results fit very well with the physical picture, 
discussed above, of the band gap closing when the cy- 
clotron radius becomes sufficiently small compared to the 
neck width of the GAL. It is also interesting to note the 
tendencies for the eigenstates to localize at the edge of 
the antidots for smaller magnetic fields, a result that fits 
well with previous studies of the magnetic properties of 
a single, isolated antidot in grapheneJ^ 



V. SUMMARY AND DISCUSSION 

Using a method based on the recursive Green's func- 
tion formalism extended to deal with structures periodic 
in two dimensions, we have calculated Hofstadter but- 
terflies of graphene antidot lattices (GALs). While the 
low-energy properties of GALs are usually well described 
in a simpler gapped graphene model, we find qualitative 
differences in the case of GALs in magnetic fields. In par- 
ticular, the lowest Landau level is not - as is the case in 
gapped graphene models - independent of magnetic field. 
Instead, we find that the GAL band gap can be effec- 
tively tuned by applying a magnetic field. In particular, 
the band gap is quenched entirely when the cyclotron ra- 
dius becomes of the order of the neck width of the GAL. 
While for lower magnetic fields, the eigenstates nearest 
the Dirac point energy are localized predominantly at 
the antidot edges, we find that in the transition region, 
where the band gap is quenched, the states become in- 
creasingly localized between antidots. Using a simple 
scaling law we show that the predicted band gap quench- 
ing might be seen for reasonable magnetic field strengths 
in experimentally feasible structures. 

We note that a similar effect of magnetically in- 
duced band gap quenching has been seen in graphene 
nanoribbonsi^ In the case of ribbons, however, the 
quenching occurs in a regime where the cyclotron radius 
is much larger than the nanoribbon width and, thus, be- 



fore the formation of Landau levels. For GALs, we ob- 
serve something quite different, namely a quenching of 
the gap that occurs for magnetic field strengths signifi- 
cantly larger than the onset of a Landau level structure. 
For ribbons, this effect has been used to predict an intrin- 
sic magnetoresistive effect;^ It is interesting to wonder 
whether very large negative magnetoresistancc might be 
seen for GALs in magnetic fields due to the magneti- 
cally induced band gap quenching. Indeed, Giesbers et 
al. have seen indications of negative magnetoresistancc 
in measurements of GAL samples with features of the 
order 100 nm in magnetic fields in the range of a few 
Tesla4 Quite large negative magnetoresistancc has also 
been seen at low temperature measurements of GALs 
with neck widths of roughly 50 nm in recently published 
results by Zhang et alM- Our scaling result would sug- 
gest a quenching of the band gap at roughly 2.6 T. which 
is of the same order of magnetic field that they see a 
significant increase in conductance. We note that we ex- 
pect disorder to play a significant role in actual samples, 
the effect of which we have not considered in the present 
paper. However, an interesting point in relation to the 
results of Fig. [7] is that the tendency of the states to 
localize between the antidots might have the benefit of 
reducing any effects pertaining to the particular edge ge- 
ometry of the hole. This would conceivably include any 
effects due to disorder on the edge of the antidot, which 
is to be expected in most realizations of GALs. We plan 
to pursuit further theoretical studies of these topics in 
future work. 
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Appendix A: Recursive method 

Our numerical method is based on standard recur- 
sive Green's function techniques, wherein advantage is 
taken of the block-diagonal form of the Hamiltonian. 
Here, we briefly outline the standard formalism and ex- 
tend it to incorporate periodic boundary conditions in 
both directions. We consider the whole system as be- 
ing assembled via N subsystems, as illustrated in Fig. |S1 
Each isolated subsystem is characterized by a Hamilto- 
nian matrix H„ and a corresponding Green's function 
G„ = \{E -\- i"f)1. — H„]~^, with I an identity matrix, 
and where we have included the broadening term 7. The 
coupling between the n'th and the n + I'th subsystem 
is described by the matrix V„. Further, we define H'"' 
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FIG. 8: The Hamiltonian H''^' of the complete system can be 
assembled from Hamiltonians of isolated subsystems Hn and 
the matrices V„ describing the coupling between the subsys- 
tems. Note that due to periodicity, the first and final subsys- 
tems are coupled via the matrix Vi_jv. 



and G^") as the Hamiltonian and Green's function, re- 
spectively, of the combined system of the first n sub- 
systems, including the coupling between them. Taking 
advantage of the block-diagonal form of H^") for n < N, 
the standard recursive Green's function method relies on 
a recursive procedure for n -I- 1 < A'^ reading 



ri(n+i) 

,("+!) 



in) 



G 



(n+l) 



(Al) 

where we have introduced the self-energy = 
Vj'jGi"^V„, which takes care of the coupling between 
the subsystems. Note that the subscripts of the Green's 
functions refer not to individual elements but to subma- 
trices defined by the sizes of the subsystems. For the last 
step of the recursive procedure, n -|- 1 = A'^, we do take 
into account the coupling to the first subsystem due to 
periodicity, resulting in 

g'^nm = HE + *7)I - Hat - l^N-iV' 



,(n+l) 



G 



(N) 
.N 



,{N-1 



a recursive algorithm specifically for calculating the den- 
sity of states. We calculate the density of states via the 
relation p{E) = — Tr^^ImlTg^''}, where we have defined 

r,^;"' = Tr[G^")]. We first treat the case where n + l < N 
or periodicity is ignored. Using the relations for the el- 
ements of the Green's function given in Eqs. (jAip . and 
doing a bit of algebra, we arrive at 

4"+^' = + Tr [g„+i f VtE(;^) V„ + i) 1 , (A4) 



where to ease notation we have defined the shorthand 
G„+i — g|j"|]^^,j_|_]^. Furthermore, we have introduced 
the ancillary matrix 



e(") _ Vg^"^g^"^ 



(A5) 



with the recursive relation 

Ef^'' = G„+i (Vt eWv„ + l) G„+i. (A6) 

While our notation is slightly different, we note that this 
result has previously been derived by MacKinnoni^i^^ 
We include the results here for completeness. For the 
final step of the recursive algorithm we take into account 
periodicity of the structure via Eqs. (jA2p which leads to 



(JV) 



G 



T, 



(N-1) 



Tr 



G 



N 



Vi, 



TV 



+ ^l,Ar-'^ln ^ N-l + N-l'^nl ^ l,N 

+V]^_,E(r'^V^_i+l)' 



(A7) 



where we have introduced additional ancillary matrices 



E. 



(n) 
nl 



E 



(n) 



E 



(n) 
11 



E„(n)„(n) 

i 
n 

Ep,(n)p,(n) 

i 
n 

EGi:-^G(:;\ 



(A8) 
(A9) 
(AlO) 



^N,j - ^N,N [ ^i.jv'^ij + ^ jv-i'-'w-ij I I J with recursive relations 



Gi,j - Gj + (^(^j 1 Vi,Ar + (j.jv-l ^N-1 j G^_^ 

(A2) 

with the self-energy 

^N-l - l,N^l,N-l + ^ N-l^N-l,N-l) ^ N-l 

+ (vt_,G(^-)+v],„,G(^-l V,M.A3) 



E 



("+!)_ 



G„ 



^nl — '^n+1 

and similarly 



Vt E(:;)+(vtE(,^,)V„+l) Gi';+i!ll ,(A11) 



E 



(«+!)_ 



While the method above can be used directly to ob- 
tain the full Green's function G'^' of the whole system, 
when dealing with very large systems this method is both 
too slow and too memory intensive. Instead, we derive 



Ei^^V^.+Gtt+l (vtEWv„+l)]G„+i.(A12) 
Finally, we have 

E(n+i) ^ El^) + Gt++^lvt eJ:;) + E(:;V„Giti!l 

+G^C+l (viEWv„ + l) Git,^}, (A13) 

(n) (n) 

while the recursive procedures for G\ „ and G„ { are of 
course given by Eqs. (jAl[) . 
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